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Q ■ We consider the neutralino as a dominant Dark Matter particle 

in the galactic halo and investigate some general issues of direct DM 
searches via elastic neutralino-nucleus scattering. On the basis of 
conventional assumptions about the nuclear and nucleon structure 
j> ; we analyse constraints on SUSY model parameter space accessible 

by the direct DM searches. This analysis shows that DM detectors 
^ . fall into the three different categories with respect to their sensitiv- 

ity to different groups of the SUSY model parameters. 

We calculate the event rate for various experimentally interest- 
Q\ ' ing isotopes within the Minimal Supersymmetric Standard Model 

(MSSM) taking into account the known accelerator and cosmolog- 
P< ical constraints on the MSSM parameter space. 

Ph! We investigate the role of nuclear spin in elastic neutralino- 

J|h | nucleus scattering. It is found that the contribution of the spin- 

dependent interaction to this process is subdominant for nuclei with 
atomic weights A > 50. 

X 

.«?. 

1. Introduction 

Analysis of the data on distribution and motion of astronomical objects within our 
galaxy and far beyond indicates presence of a large amount of non-luminous dark matter 
(for review seeS). According to estimations, dark matter (DM) may constitute more 
than 90% of the total mass of the universe if a mass density p of the universe is assumed 
to be close to the critical value p c = 3H 2 /StiGn (H is the Hubble constant and Gn is 
the Newtonian gravitational constant). The exact equality Q — p/ p c — 1, corresponding 
to a flat universe, is supported by naturalness arguments and by inflation scenarios. 
The theory of primordial nucleosynthesis restricts the amount of baryonic matter in 
the universe to ~10%. Thus a dominant component of DM is non-baryonic. The recent 
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data by the COBE satellitea on anisotropy in the cosmic background radiation and the 
theory of the formation of large scale structures of the universe lead to the conclusion 
that non-baryonic DM itself consists of a dominant (70%) "cold" DM (CDM) and a 
smaller (30%) "hot" DM (HDM) component!!. 

A possible HDM candidate is the massive neutrino. The SUSY model neutralino 
(x) is currently a favorable candidate for CDM. This is a Majorana (x c = x) particle 
with spin 1/2 predicted by supersymmetric (SUSY) models. 

There are four neutralinos in the minimal supersymmetric extension of the stan- 
dard model (MSSM) (seeu). They are a mixture of gauginos (W3, B) and Higgsinos 
(#1,2) being SUSY partners of gauge (W3, B) and Higgs (#1,2) bosons. The DM neu- 
tralino x is assumed to be the lightest supersymmetric particle (LSP) and therefore is 
stable in SUSY models with i?-parity conservation. 

In our galaxy most of the mass is expected to be in form of a spherical dark halo. 
Microlensing searches have discovered far more MACHOs in the disk than in the halo 
of the Galaxy. The data obtained from these searches is consistent with the fraction of 
MACHOs in the halo dark matter! less than 30%. Thus, most of the halo of the galaxy 
must be non baryonic cold dark matter. 

In the galactic halo neutralinos are assumed to be Maxwellian distributed in 
velocities with a mean velocity in the earth frame v ~ 320 km/se J. Their mass density 
in the Solar system is expected to be about p ~ 0.3 GeV ■ cm -3 . Therefore, at the earth 
surface neutralinos might produce a substantial flux ($ = p ■ v/M) of $ > 10 7 cm -2 
sec -1 for a particle mass of M~ 1 GeV. In view of this one may hope to detect DM 
particles directly, for instance through elastic scattering from nuclei inside terrestrial 
detector. 

The problem of direct detection of the DM neutralino x v i & elastic scattering off 
nuclei has attracted considerable efforts during the last decade and remains a field of 
great experimental and theoretical activityuiij. 

In this report we address the questions concerning prospects for the direct detec- 
tion of the supersymmetric Dark Matter with the current and the near future detectors. 

In Sec. 2 we consider general issues of such experiments for exploration of the 
SUSY model parameter space and classify possible DM detectors with respect to their 
sensitivity to different domains of this parameter space. We propose and discuss special 
criterion for assessing an isotope as a target material for a DM detector. In discussing 
general expectations for DM detection experiments, we avoid the use of specific nu- 
clear and nucleon structure models, but rather base our consideration on the known 
experimental data about nuclei and nucleon properties. 

In Sec. 3 we discuss predictions for the DM detection event rate obtained in the 
framework of the MSSM. We undertake a systematic exploration of a broad domain of 
the MSSM parameter space restricted by the well known accelerator constraints and 
by the cosmological bounds on neutralino relic abundance in the universe. The effect 
of a non-zero threshold energy of a realistic DM detector is analyzed. 

Sec. 4 is devoted to the role of nuclear spin in the DM neutralino detection. In 



2 



general, the event rate R for elastic x-nucleus scattering contains contributions from 
the spin-dependent (R s d) and spin-independent (R S i) neutralino-nucleus interactions: 
R = R s d + R S i. We have found that the R S i contribution dominates in the total event 
rate R for nuclei with atomic weight A > 50 in the region of the MSSM parameter 
space where R = R sd + R s i > 0.01 f 6 ^ - The lower bound 0.01 f t,e " fs seems to be far 
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below the sensitivity of realistic present and near future DM detectors. Therefore one 
can ignore the region where R < 0.01 as invisible for these detectors. 

In view of this result we do not expect crucial dependence of the DM event rate 
on the nuclear spin for detectors with target nuclei having an atomic weight larger than 
50. In other words, we expect essentially equal chances for J = and J ^ detectors 
to discover DM events. 

In particular, this conclusion supports the idea that presently operating (3(3- 
detectors with spinless nuclear target material have good prospects for DM neutralino 
search. 

Sec. 5 gives a conclusion. 

2. General Properties of the Neutralino-Nucleus Interactions 

A DM event is elastic scattering of a DM neutralino from a target nucleus produc- 
ing a nuclear recoil which can be detected by a detector. The corresponding event rate 
depends on the distribution of the DM neutralinos in the solar vicinity and the cross 
section a e i(xA) of neutralino-nucleus elastic scattering. In order to calculate a e i(xA) 
one should specify neutralino-quark interactions. The relevant low-energy effective La- 
grangian can be written in a general form as 



L eff = Yl ( A q ' • <?7 M 75<? + tt 1 • C q ■ XX ' QQ + O — , (1) 

q \ M W / \ m q/ 

where terms with the vector and pseudoscalar quark currents are omitted being negligi- 
ble in the case of the non-relativistic DM neutralino with typical velocities v x ~ 10 _3 c. 

In the Lagrangian (JJ) we also neglect terms which appear in supersymmetric 
models at the order of l/m| and higher, where mg is the mass of the scalar superpartner 
q of the quark q. These terms, as recently pointed out by Drees and NojiriEl, are 
potentially important in the spin-independent neutralino-nucleon scattering, especially 
in domains of the MSSM parameter space where rriq is close to the neutralino mass M x . 
Below we adopt the approximate treatment of these terms proposed in Ref.0 which 
allows "effectively" absorbing them into the coefficients C q in a wide region of the SUSY 
model parameter space. 

The coefficients A q , C q depend on the specific SUSY model and will be considered 
in the next section. 

Here we survey general properties of neutralino-nucleus {x~A) scattering following 
from the Lagrangian ([]]). 

A general representation of the differential cross section of neutralino-nucleus 
scattering can be given in terms of three spin-dependent Ti{q 2 ) and one spin-independent 
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jF 5 (g 2 ) form factors as followsH 

^(v, g 2 ) = 8 -^f{al + ooo! ■ + 4 ■ ?W) +ct-A 2 T 2 s (q 2 )) . (2) 

The last term corresponding to the spin-independent scalar interaction gains coherent 
enhancement A 2 (A is the atomic weight of the nucleus in the reaction). Coefficients 
cii, Co do not depend on nuclear structure and relate to the parameters A q ,C q of the ef- 
fective Lagrangian (JJ) and to the parameters Aq, f s , f characterizing nucleon structure. 
One has the relationships 

Go = {A u + A d ) (Au + Ad) + 2AsA s , a x = (A u - A d ) (Au - Ad) , (3) 
co = f mXu t mdCd + fsC a + hl-f s - f)(C c + C b + C t ) 

Here Aq p ^ are the fractions of the proton (neutron) spin carried by the quark q. The 
standard definition is 

< p(n)\qYlsq\p(n) >= 2S^ n) A^ n \ (4) 

where = (0, o p (n)) is the 4-spin of the nucleon. The parameters Aq p ^ can be 

extracted from data on polarized nucleon structure functioning and hyperon semilep- 
tonic decay dataS. It has been recently recognizedill that the new preliminary SMC 
measurement sEl of the spin structure function of the proton at Q 2 = 10.3 GeV 2 may 
have dramatic implications for calculations of the spin-dependent neutralino-nucleus 
scattering cross section. The values of Aq extracted from these new data in comparison 
with previous EMC@ data are much closer to SU(3) naive quark model (NQM) pre- 
dictionsi'S. This gives rise to small enhancement of the spin-dependent cross section 
for nuclei with an unpaired proton and a strong (by a factor of about 30) suppression 
for nuclei with an unpaired neutron. In view of this we use in the analysis Aq values 
extracted both from the EMC0 and from SMcEl data. 

The other nuclear structure parameters f s and / in formula (0) are defined as 
follows: 

< p(n)\(m u + m d )(uu + dd)\p(n) > = 2/M p(n) M, (5) 
< p{n)\m s ss\p{n) > = / s M p(n) M. 

The values extracted from the data under certain theoretical assumptions are§!: 

/ = 0.05 and f s = 0.14. (6) 

The strange quark contribution f s is known to be uncertain to about a factor of 2. 
Therefore we take its values in the analysis within the interval 0.07 < f s < O.30& 
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The nuclear structure comes into play via the form factors Tij{q 2 ) 1 J-'s{. < i 2 ) i n Eq. 
(0). The spin- independent form factor J-'siQ 2 ) can De represented as the normalized 
Fourier transform of a spherical nuclear ground state density distribution p(r) 



F s {q 2 ) = J ^ 3 rp(r)e irq . (7) 

In the analysis we use the standard Woods-Saxon inspired distributionil. It leads to 
the forma factor 

Ml 2 ) = 3^fir^) 2 , (8) 
qK 

where Ro = (R 2 — 5s 2 ) 1//2 and s ~ 1 fm are the radius and the thickness of a spherical 
nuclear surface respectively, ji is the spherical Bessel function of index 1. 

Spin-dependent form factors ^^(q 2 ) are much more nuclear model dependent 
quantities. The last few years have seen a noticeable progress in detailed nuclear model 
calculations of these form factors. For many nuclei of interest in DM search they have 
been calculated within the conventional shell model!! and within an approach based on 
the theory of finite Fermi systems^. We use the simple minded parametrization for q 2 
dependence of ^(q 2 ) in the form of a Gaussian with the rms spin radius of the nucleus 
calculated in the harmonic well potential^. For our purposes this semi-empirical scheme 
is sufficient. 

An experimentally observable quantity is the differential event rate per unit mass 
of the target material. It reads 



dR 



dE. 



r 





I 




J V 



dvf(v)v^(v,E r ) (9) 



Here f(v) is the velocity distribution of neutralinos in the earth frame which is usu- 
ally assumed to be the Maxwellian distribution in the galactic frame. v max = v esc ~ 
QOOkm/s is the escape velocity at the sun position; v min = (MAE r /2M 2 ed ) 1 ' 2 with Ma 
and M red being the mass of nucleus A and the reduced mass of the neutralino-nucleus 
system respectively. Note that q 2 = 2MaE t . 

The differential event rate is the most appropriate for comparing with the observ- 
able recoil spectrum and allows one to take properly into account spectral characteris- 
tics of a specific detector and to separate the background. However, for a more general 
theoretical discussion the event rate integrated over some domain of recoil energy is 
more useful and commonly employed for estimating the prospects for dark matter de- 
tection, ignoring experimental complications which may occur on the way. Moreover, 
the integrated event rate is basically less sensitive to details of nuclear structure then 
the differential one (0). The q 2 shape of the form factors J r ij(q 2 ), Fsiq 2 ) in Eq. (g) may 
essentially change from one nuclear model to another. Integration over some q 2 region 
reduces these variations drastically. 

Define the integral event rate as: 

r E 2 dR 

R(E U E 2 ) = / -—e(E max - E r )dE r (10) 
Je 1 dE r 
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Here E max = 2M^. ed v1 sc / Ma It is a common practice in theoretical papers to analyze 
the total event rate R = R(0, oo). For the realistic DM detector one should take into 
account a non-zero threshold energy E\ > E t f ir . In what follows we use the total event 
rate R in a general discussion. The effect of a DM detector threshold E t h r is analyzed 
later on in Sec. 2. 

Now let us address the question of possible constraints on SUSY models reach- 
able in the direct DM search experiments. As seen from Eq. (|2|) there are just three 
parameters ao,i; c o accumulating all SUSY model dependence via parameters A q ,C q . 
Therefore, in the experiments discussed here the only three constraints on the combined 
SUSY model parameters ao,i,co are accessible. 

For getting a more transparent physical meaning of these constrains it is useful 
to adopt the approximation of the odd group model (OGM)0. It assumes for odd-even 
nuclei that a dominant contribution to the nuclear spin comes from the odd nucleon 
group. For the most of nuclei of interest for DM searches this approximation is fairly 
good. If so, one can write the following formula for the total event rate0 



R 



■p(n) 



2 events 



W)-S(n) + ^0-C j kg day - (11) 



for the case of proton (neutron) odd group nucleus. Here, a p ( n ) = a ±ai. The parameters 
4>i depend on properties of the target nucleus as well as on the mass density and the 
average velocity of DM particles in the solar vicinity. 

The quantities a Pi „ and Co contain all dependence on the parameters of the effec- 
tive Lagrangian (|l|) and do not depend on nuclear properties. This factorization leads 
us to the following conclusions. 

From Eq . (|TT|) we see again that measuring the event rate R we can study just 
three special combinations of fundamental parameters a p n and Cq. This is the only 
information about fundamental parameters accessible in DM search experiments. R is a 
linear combination of the quantities a p n and Cq. To extract experimental limitations for 
each of them one should search for DM with different target nuclei. We can distinguish 
three categories of DM detectors with respect to their sensitivity to a Pi „ and cq. These 
are detectors built of spin-non-zero target nuclei with an odd proton (neutron) group 
probing a linear combination a p ( n ) and Cq and spin-zero target nuclei sensitive only to 
the scalar part Cq of the neutralino-nucleus interaction. 

We would like to stress again the following. To extract all possible information 
about SUSY-model parameters from direct DM search one should have three above 
mentioned types of DM detectors. No other information can be obtained from the 
direct DM search experiments. Different detectors can only improve the data on the 
three above-defined groups of SUSY-model parameters. 

Accessible constraints can be represented as experimental constraints to the ef- 
fective parameters a p n and cq accumulating all SUSY model dependence of the event 
rate. According to formula ( p|) the corresponding exclusion plots in the a p — Co and 
a n — c o planes have a form of ellipses as presented in Fig.l for nuclei with neutron 
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FIG. 1. Exclusion plots for the effective SUSY parameters a p ( n ) (Axial), 
Co (Coherent) for nuclei with odd neutron (left panel) and proton (right 
panel) groups respectively. 



and proton odd groups respectively. The plots correspond to a DM detector sensitivity 
R > 1 f ve1 ] 1 s . It is easy to see that the same picture holds for any sensitivity. The only 

kgday J * j j j 

effect is the rescaling of the a p ( n ) — Co axes. These plots allow one to assess which iso- 
tope gives more restrictive constraints at the same detector sensitivity R > Rdet- For 
instance, it is seen that among p-odd nuclei 205 T1 has potentially better prospects as 
a target material for DM detectors than 23 Na, 35 C1, 69 Ga, 127 I. Of course, one should 
remember that such a criterion is very superficial because it does not take into account 
limitations imposed by specific conditions under which an isotope can be used in a 
realistic DM detector. 

3. SUSY-model Predictions 

In order to obtain quantitative predictions for the DM detection event rate one 
should calculate the parameters A q , and C q of the effective Lagrangian ([!]) in the specific 
SUSY model. We follow the MSSM with the GUT unification conditions for gauge 
coupling constants and for soft SUSY breaking parameters. This model is specified 
by the standard SU(3) x SU(2) x £7(1) gauge couplings as well as by the low-energy 
superpotential and soft SUSY breaking terms. 

The effective low-energy superpotential is: 

W= Yl (hLHiLE + h D HiQD - huHzQfr) - 11H1H2. (12) 

generations 

Here L, E are lepton doublets and singlets; Q are quark doublets, U, D are up and 
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down quark singlets; Hi and H 2 are the Higgs doublets with a weak hypercharge 
Y = —1, +1, respectively. 

The effect of "soft" supersymmetry breaking can be parametrized at the Fermi 
scale part of the scalar potential: 

V soft = £ m\\&f + h L A L Hj,E + K B A D R 1 Qb-h v A lJ H 2 Qij- (13) 

i=scalars 

-(fiBH 1 H 2 + h.c.) 



and a "soft" gaugino mass term 

1 



Cfm — — ^ 



M X BB + M 2 W k W k + M 3 g a g a - h.c. (14) 



As usual, Mi )2) 3 are the masses of the 577(3) x SU{2) x £7(1) gauginos g, VT, B. = 
{m L , m E , itiq, mu, m D , m Hl ,m H2 } are the mass parameters of scalar fields. 

To reduce the number of free parameters we use the following unification condi- 
tions at the GUT scale M x : 

Au = A D = A L = Ao, (15) 

m L = m E = m Q = mu = m D = m Q , (16) 
Mi = M 2 = M 3 = mi/2) (17) 

9l (M x ) = g 2 {M x ) = g 3 {M x ) = g GUT , (18) 

where g 3 ,g 2 ,gi are the SU(3) x SU(2) x U(l) gauge coupling constants equal to g^ur 
at the unification scale M x . 

At the Fermi scale Q ~ M\y these parameters can be evaluated on the basis of 
the renormalization group equations (RGE)il. Equation flTTD implies at Q ~ Mw 

5 
3 



Mi = - tan 2 W ■ M 2 , M 2 = 0.3mg. (19) 



Here rrig = M 3 is the gluino mass. One can see from (]T5|)-(|T8|) that we do not exploit 
the complete set of GUT unification conditions for the soft supersymmetry breaking 
parameters, which leads to the supergravity scenario with radiative electroweak gauge 
symmetry breaking. Specifically, we do not unify Higgs soft masses m Hl , m# 2 with the 
others in Eq. (0). Otherwise strong correlations in the supersymmetric particle spec- 
trum would emerge, essentially attaching the analysis to a particular supersymmetric 
scenario. 

We analyze the Higgs sector of the MSSM at the 1-loop levefl with taking into 
account ti — tp>, — mixing between the third-generation squarks. Diagonalization 
of the Higgs mass matrix leads to three neutral mass-eigenstates: two CP-even states, 
h, H, with the masses m^, mu and the relevant mixing angle an as well as one CP-odd 
state A with the mass m^- We take the mass m^ as an independent free parameter. 
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A complete list of the essential free parameters we use in the analysis is 

tan/3,A , n, M 2 ,m A ,m ,m t . (20) 

The angle j3 is defined by the vacuum expectation values of the neutral components of 
the Higgs fields: tan/5 =< H$ > / < >. 

We fix further for definiteness m t = 174 GeV which corresponds to the CDF 
value of the top quark mass0. 

There are four neutralino mass eigenstates Xi i n the MSSM 

Xl = AfaW 3 + N a B + AT t3 H° 2 + (21) 

which are linear combinations of zero charged gauginos (W 3 , B and Higgsinos H 2 , H®). 
The unitary matrix M rotates the neutralino 4x4 mass matrix to the diagonal form. 
As usual, we denote the lightest neutralino Xi as X- 

Having specified the model one can derive the effective Lagrangian L e ff of low- 
energy neutralino-quark interactions. In the MSSM the first term of L e ff in Eq. ([!]) is 
induced by the Z-boson and q exchange whereas the second one is due to the Higgs 
boson and q exchange. 

Our formulas for the coefficients A q and C q of the effective Lagrangian take into 
account squark mixing qi — qb and the contribution of both CP-even Higgs bosons 
h, H. As pointed out in RefBil, the contribution of the heavier Higgs boson H can 
be important in certain cases. At some values of the angles ctH,/3 and the neutralino 
composition coefficients A/13, A/14 the contribution of the heavier Higgs boson H to 
the coefficients C q can be larger than the contribution of the lightest Higgs boson h. 
The above formulas coincide with the relevant formulas in Ref.li3 neglecting the terms 
~ l/m|and higher. As stated in Sec. 2, we adopt the approximate treatment proposed 
111 RefB. It allows one to take into account these terms "effectively" by introducing an 
"effective" stop quark t propagator. In the limit 9 g — > 0, where 6 q is the qi — q^ mixing 
angle, our formulas agree with Ref.B except for the relative sign between the Z and q 
exchange terms in the coefficients A q and up to the overall sign in the coefficients C q . 
These errors in Ref.@ were also observed in Ref.ll!. Now we are ready to calculate the 
event rate of neutralino nucleus scattering. 

In our numerical analysis we scan the MSSM parameter space within a broad 
domain 

20 GeV <M 2 < 1 TeV, \p\ < 1 TeV, (22) 

1 < tan/5 < 50, \A \ < 1 TeV, 

< m < 1 TeV, 50 GeV < m A < 1 TeV 

In the region where tan/? > 35 the top Yukawa dominance approximation is not 
applicable in the RGE. Therefore, we use procedure developed in Ref.il which takes 
into account the bottom and tau Yukawa couplings as well. 
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FIG. 2. Scatter plots of the total event rate R vs the neutralino (LSP) 
mass M x (left panel) and the threshold correction factor v vs the thresh- 
old recoil energy E t h r (right panel). Three curves correspond to M x = 
30, 100, 200 GeV from the bottom to the top. All plots for 73 Ge. 



Further limitations on the parameters space are imposed by the known experi- 
mental lower bounds on supersymmetric particle and Higgs boson masses from LEP 
and Tevatron measurements. 

The neutralino relic density Q x is also under control in our analysis. We calculate 
it following the standard procedure on the basis of the approximate 




!W = 2.13X1CT"U? -Nr I " -• „ ■ (23) 



Here T 7 is the present day photon temperature, T x /T 7 is the reheating factor, xp = 
Tp/M x « 1/20, Tp is the neutralino freeze-out temperature and Np is the total number 
of relativistic degrees of freedom at Tp. The coefficients a, b are determined from the 
expansion 

< a ann v >K,a + bx (24) 

of the thermally-averaged cross section < a ann v > of neutralino annihilation. We use 
an approximate treatment ignoring complications which occur when expansion ( p4| ) 
fails22l. We take into account all possible channels of the x ~ X annihilation. 

It is well known that cosmologically acceptable neutralinos should produce a relic 
density in the interval 

0.025 < Vt x h 2 Q < 1. (25) 

In this case neutralinos do not overdose the universe and account for a significant 
fraction of the halo DM. 
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As an example of our event rate calculations!!! we show in Fig. 2 (left) the total 
event rate R for 73 Ge. The scatter plot was obtained by random point generation 
in the MSSM parameter space with the constraints discussed above. More detailed 
presentation of these results including other isotopes is given in Ref.fi 

A non-zero threshold energy E r for the case of a realistic DM detector may 
essentially modify pictures displayed in Fig. 2 (left). To quantify the effect of a non-zero 
detector threshold we introduce the ratio 

RjEthr, oo) 

U{Ethr) - R(0,oo) (26) 

In Fig.2(right) this threshold factor is plotted for the isotope 73 Ge. It is seen that for the 
lighter neutralino v{E t hr) falls faster than for the heavier one. One can easily understand 
this dependence noticing that a sizable R can be obtained for E t ^ r < 10~ 6 M X . The latter 
corresponds to the mean kinetic energy of the DM neutralino. 

Combining plots Fig.2(left) and Fig.2(right) we can estimate the maximal value 
of the event rate i?(H-M) integrated from E t h r = 48 KeV which corresponds to the 
threshold recoil energy of the Heidelberg-Moscow germanium detectorE3. As seen for 
M x « 200 GeV, it approaches R(H - M) = 0.2§^ comparable with the sensitivity 
quoted by this collaboration^ Therefore, one may expect this experiment to provide in 
the near future the DM constraints on the MSSM parameter space. Detailed analysis^ 
leads us to the conclusion that in the similar position is another germanium experi- 
ment by the Caltech-PSI-Neuchatel collaboration^. Other DM experiments^ are more 
distant from the reach of the allowed part of the MSSM parameter space as yet. 

Note that in the analysis we ignore possible rescaling of the local neutralino 
density p which may occur in the region of the MSSM parameter space where flh 2 < 
0.050. This effect, if it took place, could modify the predictions for the event rate 
We assume that neutralinos constitute a dominant component of the DM halo of our 
galaxy with a density p = 0.3 GeV-cm~ 3 in the solar vicinity. 

4. The role of target nucleus spin 

To study the role of nuclear spin in elastic x- nuc l eus scattering we consider the 

ratio 

rj = R sd /R sl (27) 

characterizing the relative contribution of spin-dependent and spin-independent inter- 
actions. Here R sc i and R si are the spin-dependent and spin-independent parts of the 
total event rate R respectively. The quantity rj + 1 determines the expected relative 
sensitivity of DM detectors with spin-non-zero (J ^ 0) to those with spin-zero (J = 0) 
nuclei as target material, if their atomic masses are close in value. If 77 < 1, then 
detectors with spin-non-zero and spin-zero target materials have approximately equal 
sensitivities to the DM signal, otherwise if 77 > 1, the spin-non-zero detectors are more 
sensitive than the spin-zero ones. 
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FIG. 3. The nuclear factor t}a vs the atomic weight A 
for nuclei with non-zero spin. 



Let us consider separately the dependence of r\ on the nuclear structure param- 
eters and on the parameters of neutralino-quark interactions determined in a specific 
SUSY model. Within our approximations one can writeS 

V = TiATiSSn (28) 

The factorization of the nuclear structure i]a from the supersymmetric part of 
the neutralino-nucleus interaction r] susy is essentially based on the assumption of the 
odd group modelli3 about a negligible contribution of the even nucleon group to the 
total nuclear spin. t]a is a factor depending on the properties of the nucleus A, while 
Vsusy IS defined by the SUSY-model which specifies the neutralino composition and 
the interactions with matter. The SUSY-factor also depends on the nucleon matrix 
element parameters (^), (BD and on the shell-model class to which nucleus A belongs, 
being rjs USY for the shell-model "neutron" ( 3 He, 29 Si, 73 Ge, 129,131 Xe,...) and rf SUSY for 
the shell-model "proton" ( 19 F, 23 Na, 35 C1, 127 I, 205 T1,...). 

Fig. 3 shows the nuclear factor 77^4 versus the atomic weight ^lll. The height of 
the symbols in the picture represents the variation of the ratio tja within the explored 
interval of the neutralino mass of 20 GeV < M x < 500 GeV. 

It follows from Fig. 3 that t]a < 1 for A > 50. Thus at A > 50 there is no 
nuclear structure enhancement of the spin-dependent event rate as compared to the 
spin- independent one. 
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FIG. 4. Scatter plots of the total event rate R vs the ratio r) susy . Two 
representative nuclei with an unpaired proton ( 71 Ga) and an unpaired 
neutron ( 73 Ge) are presented. 



The next step is calculation of the SUSY-factor rfj^ y within the MSSM. We 
have performed numerical analysis of the MSSM parameter space as described in the 
previous section. The following absolute upper bound for the SUSY-factor in Eq. ( f28| ) 
was found: 

(29) 

in the subdomain of the parameter space where the total event rate is R > 0.01 et,ent . 

1 1 kg-day 

Fig.4 shows the distribution of the points in the R—r] susy plane. Plots are given for 
two representative nuclei with an unpaired proton (p-like), 71 Ga, and with an unpaired 
neutron (n-like), 73 Ge. The nuclei are taken for convenience near the point A = 50 (see 
Fig.3). For heavier nuclei we have obtained basically the same picture and our further 
conclusions correspond to all nuclei with A > 50. One can see the above quoted ( p9| ) 
upper bound rj susy < 2 for both cases presented in Fig.4. 

Now we may combine the bound ([29|) with the values of the nuclear factor tja 
represented in Fig.3. Then we obtain the conservative estimate: 

r) = Rsd/Rsi = VaVsSj ~ 1.6 for nuclei with A > 50 (30) 
at a detector sensitivity up to R > 0.01 ^ e ^ . However, as is seen in Fig.4, the majority 



of points generated in the domain (|22| ) of the MSSM parameter space are concentrated 
at r] < 1. The tendency is that at higher sensitivities (lower R accessible) we get 
i] < 1 for heavier nuclei and vice versa. 

In Fig. 5 we also present plots of the ratio r(A) = R sd (A) / R sd (J 3 Ge) of the spin- 
dependent part R s d of the event rate for some target materials (A) to that for 73 Ge. 
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LSP Mass LSP Mass 

FIG. 5. The ratio r(A) = R sd (A)/R sd ( 73 Ge) of the spin-dependent 
event rate R sd for nuclei 19 F and Nal to the spin-dependent event rate 
R sd ( 73 Ge) for 73 Ge. 



This ratio for 129 Xe is r(_^ 9 Xe) ph 1.2 being fairly independent of the neutralino mass 
M x . As explained in Ref.cJ this is the case because both 73 Ge and 129 Xe are nuclei with 
an unpaired neutron. 

It follows from our detailed analysis^! and is illustrated by Fig. 5 that the maximal 
values of the total event rate for Nal, 73 Ge and 129 Xe are typically the same while for 
CaF2 they are lower by about a factor of 5. On the other hand, the sensitivity of CaF2 
to the spin-dependent part of the neutralino-nucleus interaction is by about a factor of 
10 larger than that of Nal, 73 Ge and 129 Xe. The last three materials have approximately 
an equal spin sensitivity. 



5. Conclusion 

The central point of this report is that the operating Dark Matter detectors are 
in the position to probe the MSSM parameter space in the near future. We argued that 
there are the only three sorts of constraints attainable in direct DM searches. In this 
respect DM detectors fall into three different categories probing different combinations 
of SUSY model parameters. They are detectors built of target nuclei with zero spin 
and with non-zero spin having odd proton, odd neutron groups. 

It was pointed out that for sufficiently heavy nuclei with atomic weights A > 50 
the spin-independent event rate R S i is typically larger than the spin-dependent one 
R sd if low rate DM signals with total event rates R = R sd + R S i < 0.01 are 
ignored. This cut-off condition reflects the realistic sensitivities of the present and the 
near-future DM detectors. 

The main practical issue is that two different DM detectors with (J = 0, Ai) and 
with (J 7^ 0, A2) nuclei as a target material have equal chances to discover DM events 
if Ai ~ A 2 > 50. 

Another aspect of the DM search is the investigation of the SUSY-model param- 
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eter space from nonobservation of DM events. For this purpose experiments both with 
J = and J ^ nuclei are important. 

We have compared several examples of popular (see for instance0 and references 
therein) materials with non-zero spin nuclei as a target in a DM detector. We have 
not found an essential difference between Nal, 73 Ge and 129 Xe as a target material 
for DM detectors from the point of view of their total and spin sensitivity. We expect 
these materials to have better prospects as compared with CaF 2 for discovering DM 
events. The former materials have a total event rate by about a factor of 5 larger than 
the latter one. On the other hand, CaF2 can give a more stringent constraint on the 
spin-dependent part of the event rate, having a spin sensitivity by about a factor of 10 
larger than Nal, 73 Ge and 129 Xe. 

Estimating prospects for various isotopes as target materials of DM detectors 
one may use the criterion given in Sec. 1 in the form of exclusion plots (Fig.l) for the 
effective SUSY model parameters. These plots characterize sensitivity of an isotope to 
the DM event ignoring possible complications of their employment in a detector. 

In conclusion we would like to stress that the efforts to improve existing and 
to construct new DM detectors are justified by the expectation that these detectors 
will be able to tackle supersymmetry from the side unreachable for the accelerator 
experiments. Therefore, these two types of experiments are complementary and should 
be considered as equally important in searching for supersymmetry. 
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